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ABSTRACT 

The black hole candidate XTE J1908+094 went into outburst for the first time since 2003 in October 
2013. We report on an observation with the Nuclear Spectroscopic Telescope Array (NuSTAR) and 
monitoring observations with Swift during the outburst. NuSTAR caught the source in the soft state: 
the spectra show a broad relativistic iron line, and the light curves reveal a ^40 ks flare with the count 
rate peaking about 40% above the non-flare level and with significant spectral variation. A model 
combining a multi-temperature thermal component, a power-law, and a reflection component with 
an iron line provides a good description of the NuSTAR spectrum. Although relativistic broadening 
of the iron line is observed, it is not possible to constrain the black hole spin with these data. The 
variability of the power-law component, which can also be modeled as a Comptonization component, 
is responsible for the flux and spectral change during the flare, suggesting that changes in the corona 
(or possibly continued jet activity) are the likely cause of the flare. 

Subject headings: accretion, accretion disks - black hole physics - stars: individual (XTE J1908-I-094) 
- X-rays: binaries 


1. INTRODUCTION 

XTE J1908-I-094 is an X-ray transient serendipi- 
tously discovered with the Rossi X-ray Timing Explorer 
(RXTE) Proportional Counter Arr ay (PCA) when i t 
went into outburst in 2002 February ([Woods et al.ll200§) . 
The source flux i n the 2--10keV band rose by a factor of 
^3 in one month ([Woods et al.ll2002h and reached about 
100 mCrab on 2Ci02 April 6 (iGogiis et al]l20(34ll . The 
power density spectrum showed a broad quasi-periodic 
oscillation (QPO) at 1 Hz witho ut any coherent pu lsa- 
tion between 0.001 and 1024 Hz ([Woods et al.l[2?)f)^ . In 
the X-ray energy spectrum, an iron emis sion line and 
a hard tail up to 250 keV were detected (iWoods et al.l 
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l2002t iFeroci et al.ll200Jl . The hard X-ray spectrum did 
not agree with an extrapolation of the abs orbed power- 
law rn odel in the 2-30 keV band reported bv IWoods et al.l 
( 200211. and showe d a high-energy cut-off at 100 keV 
( Feroci et al.ll200Jl . Given the timing and spectral char¬ 
acteristics, XTE J1908-I-094 is suggested to be a black 
hole candidate (IWoods et al.l l20(3^ IFeroci et al.l 120021 : 
lin’t Zand et al.l[2002t iCdgiis et al.ll2004[l. 

The broad band X-ray spectrum of XTE J1908-I-094 
is well htted with two continuum components, a multi¬ 
temperature disk blackbody with kT ~ 0.8 keV and a 
Compton plasma with a temperature near 40 keV, and 
an emiss ion line centered on the location of the Fe 
Ka line ([in’t Zand et al.l [2002h . The emission feature 
is very broad with FWHM= 3.2 ± 0.5 keV (line width 
cr = 1.4 ± 0.2 keV), which may be due to Compton scat¬ 
tering in a corona or the relativistic effects from gravita- 
tional redshift and Do ppler broadening of orbital motion 
([in’t Zand et al.ir2002[l . In the latter case, the broadening 
of the Fe emission line would be expected to be asymmet¬ 
ric, a nd could be used to mea s ure the blac k hole (BH) 
spin (iRevnolds fc Now^ 20031: iMilled 1200711 . Based on 


Miller et al 


( 200911 measured the dimen- 


this method, 

sionless spin of XTE J1908-I-094 to be a = 0.75 ± 0.09. 

The radio counterpart of XTE J1908-I-094 was 
discovered with the Very Large Array (VLA) 
at R.A .=19‘®08"®53."077. Decl.=-f09°23'04."90 

(J2000.0) ([Riipen et al.l l2002tl . which is consistent 


with the Chandra position, R.A. =1 9®®08™53.’^07, 
Decl^-|-09°23'05."0 ([Jonker et al.l 1200411 . iChatv et al.l 
([2002f l identified a likely near-infrared (NIR) counter¬ 
part for the source, but the possible counterpart was 
resolved into two sources separated by ^ 0.8arcsec in 
subsequent observations ([Chatv et al.l l200(ill . Both po¬ 
tential counterparts are consistent with XTE J1908-I-094 
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being a low mass X-ray binary (LMXB). One of the 
NIR sources would indicate an intermediate/late type 
(A-K) main-sequence companion star, while the other 
would indicate a late-type main-s equence companio n 
star with spectral type later than K (jChatv et al.ll200^ . 
Recently, Swift/\JVOT observed the X-ray source; 
however, no counterpart was found in the V band with 
the 3a limiting magnitude to be R >20.3 (|Krimm et al.l 
l2ni3bH . 


The distance to XTE J1908-I-094 is not well estab¬ 
lished. Based on its X-ray flux, the sou rce is suggested 
to be at a distance greater than 3kpc (lin’t Zand et al.l 
1200211 . From the optical measurements, th e possible dis¬ 
tance range is 3-10 kpc (|Chatv et al.ll2006ll . An estimate 
using the X-ray and r adio fluxes puts the sou rce at a 
distance of ^ 2-10 kpc (jMiller- Jones et ;ani2013[ l. 

Previously, XTE J1908-I-094 went through two out¬ 
bursts in 2002 and early 200 3 with very similar spectral 
evolution (|Gogus et al.ll2004l . On 2013 October 26, an¬ 
other outburst of XTE J1908-I-094 was detected by the 
Swift Burst Alert Telescope (BAT), with the 15-50keV 
flux re aching ~ 60 mCrab two days later (|Krimm et al.l 
l2013all . Subsequentl y, a number of telesc opes, includ- 
i ng NuSTAR, Swift (iKrimm et al.l I2013alf0 l. the VLA 
(|Miller-Jones et all 1201311 . t he Arcminute Microk elvin 
Imager (AMI) Large Array (|Rushton et al.ll2013all. the 
Moni tor of All-sky X-ray Image (MAXI. iNegoro et al.l 
[2?fl3ll a nd the Australia Telescope Compact Array 
(ATCA, iCoriat et al.l 120131 1 carried out observations of 
the source. In this paper, we report on the NuSTAR and 
Swift observations of the 2013 outburst in detail (Sec¬ 
tion [2]) and investigate its spectral evolution and prop¬ 
erties (Section [3]). We present a discussion of the results 
in Section |4] and conclusions in Section [S] 


2. OBSERVATIONS 

XTE J1908-I-094 was monitored with many short 
observations by th e Swift/X-iay Telescope (XRT) 
(|Burrows et al.lf2005ll from 2013 October 26 to 2 013 De¬ 
cember 3, and NuSTAR (jHarrison et al.l 1201311 carried 
out an observation with an effective exposure time of ^ 
45 ks on 2013 November 8. In order to study the outburst 
properties of the source, we used the NuSTAR observa¬ 
tion and all of the Swift/XKT observations that were 
long enough to achieve sufficient statistical quality (see 
Table [U for the observation list). 

2.1. NuSTAR 

The NuSTAR data (ObslD 80001014002) were pro¬ 
cessed using version 1.3.1 of the NuSTARDAS pipeline 
with NuSTAR CALDB version 20131223. The spectra 
and light curves were extracted from a region centered 
at the position of XTE J1908-I-094 with a radius of 120". 
The source region was contaminated by stray light from 
the nearby bright source GRS 1915-1-105. Thus, the 
background region was chosen carefully. We used a cir¬ 
cular background region with a radius of 80" from the 
part of the field of view that was illuminated by the 
GRS 1915-1-105 stray light and as far away from XTE 
J1908-I-094 as possible. The background count rate is 
less than 6% of the source count rate, which means that, 
even considering the stray light, the source still strongly 
dominates the spectra and light curves. The spectra of 


TABLE 1 

NuSTAR AND Swift Observations 


ObsID 

Observed date 

Exposure (s) 

NuSTAR 

80001014002 

2013-11-08 

45061/45404“ 

Swift XRT 

00033014001 

2013-10-29 

972 

00033014002 

2013-11-01 

1224 

00033014003 

2013-11-03 

1032 

00033014004 

2013-11-08 

1044 

00033014005 

2013-11-09 

1007 

00033014006 

2013-11-10 

503 

00033014007 

2013-11-11 

880 

00033014008 

2013-11-12 

975 

00033014009 

2013-11-13 

1343 

00033014010 

2013-11-14 

792 

00033014011 

2013-11-15 

1085 

00033014012 

2013-11-16 

1150 

00033014014 

2013-11-18 

979 

00033014015 

2013-11-19 

1008 

00033014016 

2013-11-20 

978 

00033014017 

2013-11-23 

958 

00033014018 

2013-11-28 

824 

00033014019 

2013-12-03 

976 


Note. — “ The exposure times of NuSTAR FPMA and FPMB, 
respectively. The Swift XRT data are taken in windowed mode. 


the two NuSTAR focal plane modules A and B (FPMA 
and FPMB), were rebinned to have at least 50 counts per 
bin. The light curves were binned to a time resolution of 
100s. 


2.2. Swift 

We reduced the Swift/XKT data from 2013 October 29 
to 2013 December 3 (see Table [ij. All data were taken 
in windowed timing mode. Using XSELEGT with XRT 
CALDB version 20140709, the spectra were extracted 
from a circular region with a radius of 20 pixels (~ 47"). 
The background extraction region is a box 20 pixels long, 
centered 100 pixels from the middle of the source extrac¬ 
tion region. Ancillary response files were created using 
the ftool xrtmkarf. At lower energies, the windowed 
timing mode shows a bump between 0.4-1 keV and a 
turn up at the lowest energie^. In order to reduce the 
low-energy spectral residuals, the grade 0 data and the 
position-dependent response matrice^^ from the latest 
XRT calibration files were used. Finally, the extracted 
spectra were rebinned to contain a minimum of 25 counts 
per bin. 


3. RESULTS 
3.1. Swift Monitoring 

The Swift monitoring observations reveal a clear evo¬ 
lution starting from 2013 October 25 (MJD 56590) 

(see Figure [T|). The Swift/BAT count rate in the 15- 
50keV bancO increased rapidly from 0.0022 ± 0.0008 
cts cm“^ s“^ on MJD 56590 to 0.026±0.002 cts cm“^ s“^ 
on MJD 56595 and then decreased sharply to ~ 0.0015 
cts cm“^ s“^ and stayed close to that level after MJD 
56604. In Swift/XRKs 0.3-10 keV band, the source 
brightened from 9.2±0.1 cts s“^ on MJD 56595, reaching 

http://www.swift.ac.uk/analysis/xrt/digest_cal.php 
http://www.swift.ac.uk/analysis/xrt/rmfs.php 

Available: http://swift.gsfc.nasa.gov/results/transients/weak/XTEJ1908p094 
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Fig. 1.— From top to bottom, the 15-50 keV Swift/3AT count rate, 
the 0.3-10 keV Sm/t/XRT count rate and the hardness ratio between 
the Swift/ym^T hard band (2.5-10 keV) and soft band (0.3-2.5 keV). 
Two vertical dashed lines indicate the time boundaries of the NuSTAR 
observation. 
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its peak count rate of 36.2 ± 0.2 cts s“i^ on MJD 56607 
and then dimmed. The hardness, defined as the ratio of 
the count rates in the 2.5-10 keV to 0.3-2.5keV bands, 
started to decrease from 2.15 ± 0.05 on MJD 56595 to 
1.030 ± 0.011 on MJD 56605, and then stayed at a value 
of 1. All of these measurements suggest that the source 
entered a state transition around MJD 56595 and was in 
the soft state 10 days later. The long exposure obtained 
with NuSTAR between MJD 56605 and MJD 56606 oc¬ 
curred after the source reached the soft state. 

First, we fitted the 0.5-10keV Swift/NBH: spectra us¬ 
ing a single absorbed power-law model. The Swift data 
below 0.5 keV were ignored during the spectral fits in or¬ 
der to exclude the low-energy spectral residuals in win¬ 
dowed timing mode. The values of the photon index, T, 
and the reduced plotted in Figure [2^. The value 

of r increased steeply from 1.6 on MJD 56595 to 4.4 on 
MJD 56605, and remained at ~ 4.5 until MJD 56629, 
consistent with the source going through the hard to 
soft state transition. After the source begins the state 
transition, the accretion disk is significant for most of 
the observation. For these observations, a single power- 
law does not provide a good fit to the spectra, and the 
addition of a disk-blackbody component provides a sig¬ 
nificant improvement to the fit. The inner disk temper¬ 
ature Tin, the normalization of the diskbb model, the 
photon index F, and the reduced shown in Fig¬ 

ure The absorbed disk blackbody plus power-law 
model could successfully fit all spectra, with Tin increas¬ 
ing from 0.3 before the state transition and stabilizing at 
about 0.7-0.8keV in soft state. 



MJD-56000(days) 


Fig. 2.— The evolution of 5ii;i/t/XRT spectral parameters, (a) photon 
index F and reduced when fitting with a single absorbed power-law 
model; (b) photon index F, inner disk temperature Tin, normalization 
of the diskbb model and reduced x^ when fitting with a two-component 
model consisting of power-law and disk components. The arrows indi¬ 
cate the upper limit of F where the lower error bars of F could not be 
well constrained. Two vertical dotted lines show the time interval of 
the NuSTAR observation. 


3.2. NuSTAR Spectroscopy 

The NuSTAR light curves of FPMA and FPMB (see 
the top panels of Figure El with background subtraction 
show a flare of ~ 40 ks duration with the peak rate being 
^40% above the non-flare rate. The background light 
curves are also shown in Figure El in order to evaluate 
if the variability might be from the nearby source GRS 
1915-1-105 rather than XTE J1908-I-094. The background 
light curves are stable at an average value of 1.3 cts s“i, 
less than 6% of the net source count rate. Thus, although 
the high background caused by GRS 1915-1-105 affects 
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Fig. 3.— The NuSTAR observation of XTE J1908+094. Top: the 3— 
79 keV light curves of XTE J1908+094 observed by FPMA and FPMB, 
respectively. Middle: the background light curves in the 3—79 keV 
band. Bottom: the hardness ratio defined as the ratio of the count 
rates in the 10—79 keV to 3—10 keV bands. The two vertical dash-dotted 
lines exhibit the duration of the fiare, the two vertical dashed lines show 
the time interval of Swift observation ObsID 00033014004, and the two 
vertical dotted lines indicate Swift observation ObsID 00033014005. 


the statistical quality of the XTE J1908+094 spectrum, 
Figure 13] demonstrates that the flare in the light curves 
comes from XTE J1908+094. To study whether the flare 
has a different spectrum from the non-flare emission, we 
first checked the ratios of the 10-79 keV count rates to 
the 3-10 keV count rates. During the flare, this hardness 
ratio increased (bottom panels of Figure E]), indicating 
that there is spectral variation. 

To investigate further, we extracted the 3-79 keV spec¬ 
tra prior to the flare, during the flare and after the flare 
and fitted them together using a simple model combining 
an energy-independent multiplicative factor (constant), 
an absorptio n mode l (tbabs), adopting abundances from 
iWilms et al.l (gOOO), a power-law model (pegpwrlw) and 
a multi-temperature disk-blackbody model (diskbb), i.e, 
constant * tbabs * (pegpwrlw -|- diskbb). Untying 
the model parameters individually or in combination, we 
found that only changing the power-law model could ex¬ 
plain the variability, with a reduced = 1-25 for 2697 
degrees of freedom (dof). As shown in the top panel of 
Figure m the power-law component changes significantly 
between the flare and non-flare spectra: before the flare, 
r = 1.96, during the flare, T = 2.23, and after the flare, 
r = 2.03. Here, we quote the best fit parameters without 
error bars because this simple model does not provide an 
acceptable fit to the data. Moreover, as shown in Fig¬ 
ure Sb, c and d, all spectra exhibit similar residuals when 
the power-law parameters are allowed to vary prior to the 
flare, during the flare and after the flare. Very poor hts 



Energy (keV) 

Fig. 4.— The NuSTAR spectra, model and residuals, (a) the unfolded 
NuSTAR spectra and model in * /(^) representation. The model 
includes a fixed disk blackbody component and a free power-law model 
before, during and after the fiare. Black and red crosses are NuSTAR 
FPMA and FPMB spectra before the fiare, respectively; green and 
blue crosses indicate the spectra during the flare; cyan and magenta 
crosses indicate the spectra after the fiare. The top dotted line exhibits 
the disk blackbody component, and the three lower dashed lines show 
the power-law components during the different stages. The power-law 
component is stronger during the flare, while the power-law compo¬ 
nents before and after the fiare show fluxes that are lower and similar 
to each other. Also, the power-law index of the flare is softer than those 
of the non-flare spectra, (b) - (d) Data/model ratio before, during and 
after the flare, respectively. 


are obtained if the power-law component is required to 
be the same for all three spectra. All of this suggests 
that the corona, rather than the mass accretion rate and 
the accretion disk, went through great changes during 
the NuSTAR observation. 

Given the strong spectral variability during the flare 
and the similar spectral properties before the flare and 
after the flare, the NuSTAR data in the 3-79 keV band 
was divided into two parts: the flare spectra and non¬ 
flare spectra. The two Swift/NKT observations, ObsID 
00033014004 and 00033014005, from prior to the flare 
and after the flare, respectively (see FigureE]), were com¬ 
bined with the non-flare spectra. Then, we used the 
model constant * tbabs * (pegpwrlw -|- diskbb) (model 
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TABLE 2 

Spectral fitting of XTE J1908+094: part I 


Model 

c 

JVh 

Ti 

Ts 

A’pli 

A/^PL2 

-^centl 

Tjcent2 

0-1 

0-2 

-^gaussl 

N gauss2 

kT 

N disk 

x"/dof 

1 

0.993 

1.6 

2.00 

2.23 

385 

554 




0.783 

650 

2709.8/2099 

2a 

0.993 ± 0.003 

2.5* 

1.95 ± 0.03 

2.15 ± 0.03 

387 ± 4 
545t^ 

z? oo-f-0.14 

.-HO.3 

o.98j;°;J^ 

. O + 0.3 
-^•'^-0.2 

'^-0.14 

o rj-H 1 • 0 
^•^-0.6 

0.755 ± 0.003 

8731™ 

2321.7/2094 

2b 

0.993 ± 0.003 

r O-I-0.8 

^•O_0.6 

1.99 ± 0.03 

2.16 ± 0.03 

391 ± 4 

551 ± 6 

r O-1-0.6 
^•'^-1.0 

4 Q+0-5 

1 c-l-0.3 
-^•0-0.2 

1 g+0.3 

lofl 

U.U< J-_0.026 

2300j;®““ 

2207.8/2093 


Note. — 

Model 1: constant * tbabs * (pegpwrlw + diskbb). 

Model 2: constant * tbabs * (pegpwrlw + gaussian + diskbb); 2a: fixed A^h at 2.5 X 10^^ cm~^; 2b: A^h was set as a free parameter. 
Model 1 is not an acceptable fit to the spectrum so we just quote the best fit parameters without error bars. 

*: fixed value; 

C is the NuSTAR FPMB normalization factor relative to FPMA; 

A^H is the X-ray absorption column density in units of 10^^ cm“^; 

Pi and T 2 are the power-law photon indices of the non-flare and flare spectra; 

A^pli and A^pl 2 are the power-law component flux normalizations over the 3—79 keV energy band in units of 10”^^ erg cm“^ s“^; 
T^centi and Ecent 2 are the Gaussian emission line energies in keV; 

(7i and <72 are the line widths in keV; 

A^gaussi and A^gauss 2 are the Gaussian component normalizations in units of 10~^ photons cm~^ s“^; 
kT is the accretion disk temperature of the diskbb model in units of keV; 

-^diskbb is the normalization of the diskbb model; 

All errors and limits are at 90% confidence level. 
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Fig. 5. — The residuals for the best fit in different models. Black and 
red symbols are NuSTAR FPMA and FPMB spectra in the non-flare 
state; green and blue symbols are the spectra of the flare. 


1 ) to fit the flare spectra plus combined non-flare spec¬ 
tra, and freed the power-law component in these two data 
sets. We And that Swift/'KBH: and NuSTAR have residu¬ 
als that are not consistent with each other in the soft X- 


ray region where they overlap. The residuals are also not 
the same for the two Swift observations. Note that the 
exposure times of the Swift observations are about 1 ks 
(see Tabled]), much shorter than that of NuSTAR, allow¬ 
ing for the possibility that Swift might catch short-term 
spectral variations in its short snapshots. For NuSTAR, 
the largest residuals are in the iron Ka emission line re¬ 
gion (Figure |4 ]d, c and d) rather than in the soft X-ray 
band observed by Swift. Therefore, in the following, we 
fit the NuSTAR spectra alone. 

As shown in Figure |3| and [S] a strong reflection com¬ 
ponent is apparent in the residuals of this fit (model 
1), leading to a large reduced = 2709.8 for 2099 
dof (see Table [5]). Similar to some other Galac¬ 
tic X-ray binaries, the reflection component is com¬ 
posed of an ir on Kg emission line and a broad re¬ 
flection excess (iLigFtman fc White! Il988t iMilleij l2007t 
iTomsick et all 1201411 . The emis sion line feature was 
also detected in th e 2002 out burst (lin’t Zand et alJl2002l : 


Gogfis et al. 


1200411 Following lin’t Zand et alJ ( 2002h and 


Gogfis et al.l (|2004ll . we used the Gaussian emission line 
model gaussian to fit this feature and performed fits 
with the neutral hydrogen column density. Ah, fixed to 
2.5 X 10^^ cm~^ (model 2a). We also tested fits where 
Ah was a free parameter (model 2b). Adding a Gaussian 
significantly improves the spectral fits with > 400 
(see Table [2] and Figure [5|). The unabsorbed disk flux 
fractions, i.e., the relative disk flux contribution to the 
total, unabsorbed flux in the 2-20 keV range, are larger 
than 80% for both the flare an d non-flare spectra, which 
meet the soft state criterion of iMcClintock &: Remillai^ 
(1200611 and also confirm that the NuSTAR observation 
was taken in the soft state. The measurement of Gaus¬ 
sian line centroid, Ecent, is dependent on Ah. Freezing 
Ah at 2.5 x 10^^ cm“^, Ecent is in the iron line region (6.4 
- 7.1 keV); leaving Ah as a free parameter, Ecent is well 
below this energy region. Given this, we then tested fits 
with Ah fixed at 4.3 x 10^^ cm“^, the average Ah when 
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TABLE 3 

Spectral fitting of XTE J1908+094: part II 


Para. 

Model 3a 

Model 3b 

Model 4 

Model 5 

c 

0.993 ± 0.003 

0.993 ± 0.003 

0.993 ± 0.003 

0.993 ± 0.003 

JVh 

4.1 ± 0.3 

4.4 ± 0.3 

5 1+°'^ 

^•-*--0.3 

5.3 ± 0.4 

El 

-L-O‘^_0.06 

2 qi+0-04 

^•Ci_o 05 

2.02 ± 0.04 

0 no+^'^^ 
^■^‘^-0.07 

Es 

^•CZ_0.05 

O 1 r4-0.04 

^•-LO_o.05 

2 10+0.05 

-^•-LO_0.04 

2.20 ± 0.04 

-Efoidi 

100* 

500* 

500* 

500* 

-Efold2 

100* 

500* 

500* 

500* 

-NpLl 

Q^g 

0 QQb+O Ol? 

U.UOO_o 024 

0.0361“;“)? 


-NpL2 

U.UO_o 03 

U.UO_o 03 

0.06 ± 0.06 


/scati 




0.010 ± 0.010 

/scat2 




U.UJU_Q Q^2 

kT 

0.719 ± 0.007 

0.711 ± 0.007 

0.6891?'“?? 


-Ndiskbb 


14201:)™ 

18101™? 


Mbh 




O 0 + 15.7 

M 




1 0+6.6 
-*-"^-0.7 

E>bh 




10* 

A^kerrbb 




1 7+4.9 
' -0.8 

5i 

5300t)™o 

4200l)™« 

9000'^™'^ 

yUUU_3QQQ 

57001)??? 


10400 ± 1700 

10000 ± 2000 

1 9500'^®°° 

iyouu_ggQQ 

119001)??? 

-Nrefl 

1 o o-l-0.16 
-L-AO_o.i5 

2 4+*^-^ 

^•^-0.7 

1 0 + 0.5 
-*^•^-0.6 

^•^-0.3 

A'ref2 

1.46 ± 0.14 

2 4+°® 

1.2 ± 0.5 

2.2 ± 0.3 

Fe/solar 

1.5* 



1.5* 

a 



-0.9981?'® 

-0.961?;?? 

i 



27l? 

33l? 

x’^/dof 

2256.9/2095 

2227.7/2094 

2209.5/2092 

2208.5/2092 


Note. — 

Model 3: constant * tbabs * (reflionxjic + cutoffpl + diskbb); 3a: 
fix Efoid — 100 keV and Fe/solar — 1.5; 3b: fix Efoid — 500 keV and 
thaw Fe/solar. 

Model 4: constant * tbabs * (relconv * reflionxjic -f- cutoffpl + 
diskbb). 

Model 5: constant * tbabs * (relconv * reflionxjic + simpl * kerrbb). 
Efoidi -Efoid 2 S'l'e the folding energy of exponential rolloff for the 

non-flare and flare spectra in units of keV; 

-NpLi and ArpL 2 are the cutoff power-law normalizations at 1 keV in 
photons keV~^cm~^s“^; 

/scati and /scat 2 are the scattered fractions of the simpl model ; 

Mbh is the black hole mass in units of the solar mass; 

M is the disk mass accretion rate in units of 10^® g sec“^; 

Dbu is the distance of the black hole in units of kpc; 

A^kerr is the normalization of the kerrbb model; 

^1 and ^2 are the ionization parameters of the reflionxjic model in 
units of erg cm s“^; 

Nj-efi and A/ref 2 are the normalizations of reflected spectrum 
(reflionxjic) in units of 10“^; 

Fe/solar is the abundance of iron relative to solar value; 
a is the dimensionless black hole spin; 

i is the inclination angle of the accretion disk in units of degree; 
other parameters are the same as in Table All errors and limits are 
at 90% confidence level. 


fitting the Swift spectra in soft state with a two compo¬ 
nent model consisting of power-law and disk components 
(Section 13.11) . We obtained Ecenti = 6.21 q 3 keV and 
£'cent 2 = 5.9±0.3 keV, and the line widths cti = 1.29;';Q;i4 
keV and (T 2 = 1.Simons keV, respectively, for the non- 
flare and flare spectra, with y^/dof = 2224.1/2094. 

Instead of the Gaussian emission line model, we then 
used the more physical model reflionxjic to fit the 
reflection component, and replaced the simple power- 
law model by a power-law with an exponential cut¬ 
off cutoffpl (model 3). The refl ionxjic model is 
an update of the m odel reflionx (|Ross et al.l 119991 : 
iRoss fc Fabianl[2?)?)^ . which calculates the reflected spec¬ 
trum from an optically thick atmosphere ionized by illu¬ 
minating X-rays with a cutoff power-law spectrum. The 
power-law photon index of reflionxjic is linked to that 
of cutoffpl. Compared with reflionx, the folding en¬ 
ergy HighECut in reflionxjic is a free parameter also 
linked to that of cutoffpl. In addition, the ionization 
parameter, and the abundance of iron, Fe/solar, ex¬ 
tend over larger ranges in reflionxjic. 

When left as a free parameter, the best fit value for 
the exponential folding energy, HighECut, is 500keV, 
which is the upper limit of the parameter range. As 
this parameter is not well-constrained, we performed fits 
with HighECut fixed at 100 keV and 500 keV, respec¬ 
tively. Moreover, we also performed fits with Fe/solar 
fixed at the initial value of 1.5 and as a free parameter. 
Good hts with reduced less than 1.08 were obtained if 
a reflection component was added. Changing HighECut 
from 100 keV to 500 keV, or unfreezing Fe/solar, other 
model parameters change only slightly, as seen for model 
3a (HighECut = lOOkeV, Fe/solar = 1.5) and model 3b 
(HighECut = BOOkeV, free Fe/solar) in Table[3l Similar 
to the power-law photon index, the ionization parameter 
in the flare stage is also larger than that in the non-flare 
stage. 

The iron Ka emission line may be distorted by rela- 
tivistic effects; ther efore, a convolution model, relconv 
(|Dauser et al JMg), was adopted to calculate relativis¬ 
tic smearing (model 4). The relconv model also allows 
for a broken power-law emissivity function for the inci¬ 
dent emission. Compared with other relativistic smear¬ 
ing models, relconv extends the black hole spin param¬ 
eter range to negative values, corresponding to a disk 
rotating counter to the black hole’s spin. 

The fits also favored a high folding energy and were 
performed with HighECut fixed to 100 keV and 500 keV. 
We included fits with the iron abundance free and also 
fixed to a value of 1.5 solar. Similarly to before, freezing 
Fe/solar or changing HighECut causes little difference 
in the residuals and other model parameters. The inner 
disk radius was set to be at the innermost stable circular 
orbit (ISCO), and the outer disk radius was set to 400 
rg, where Vg = GM/c^ is the gravitational radius. The 
emissivity indices were fixed at the default values, and we 
noted that thawing these parameters or hxing the inner 
emissivity index at 3 < gin < 10 and the outer emissivity 
index at 0 < 9out < 3 (e.g., gin = 5 and gout = 2, or 
gin = 8 and gout = 1) did not improve the fits signifi¬ 
cantly (the decrease in was less than 2.7). The best 
fit model is shown in Table [3] and Figure [S] Adding a 
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Fig. 6.— The unfolded NuSTAR spectra and components of model 
5 in * f{^) representation. Black and red crosses are, respectively, 
FPMA and FPMB spectra from the non-flare times; green and blue 
crosses are the flare spectra. The two upper dotted lines exhibit the 
simpl * kerrbb components of the flare and non-flare spectra, while the 
two lower dashed lines exhibit the reflection components. The reflection 
component of the flare spectra is stronger than that of the non-flare 
spectra. 


relativistic blurring model led to only a marginally sig¬ 
nificant improvement in For the spin of the black 
hole, a wide range is allowed, with the full parameter 
range (from —0.998 to 0.998) being covered when all the 
models we used are considered. This will be discussed in 
Section |4l 

In order to constrain the spin of black hole, the diskbb 
model was repl aced by a mor e physical disk blackbody 
model, kerrbb (jLi et al.ll2005l) . The model calculates the 
disk continuum around a Kerr black hole and fully takes 
the relativistic effects into account. More over, following 
previous papers (e.g., iTomsick et al.l[20Tll . an empirical 
Comp tonization convolution model, simpl (jSteiner et al.l 
120091) . which assumes that a fraction of seed photons are 
scattered into a power-law component, was used instead 
of the power-law model (model 5). 

Similar to the fits above, a high folding energy was 
preferred by model 5. Although we also tested the fits 
with Ffoid fixed at 100 keV and the iron abundance left 
as a free parameter, we only show the spectral fitting 
with Afoid = 500 keV and Fe/solar=1.5 in Tableland 
FigureObecause there is only a slight change in the good¬ 
ness of fit for other values of those parameters. The dis¬ 
tance of XTF J1908-I-094 is thought to be ^ 2 — 10 kpc; 
thus, Dbh was set to be 2 kpc or 10 kpc. The spin and 
the inclination of kerrbb are linked to those of relconv. 
Other model parameters were fixed at the default val¬ 
ues. We obtained a very small improvement in the fits 
with the reduced of 1.06 for 2092 dof Except for 
the normalization of the kerrbb model, all model pa¬ 
rameters show little changes if Dbh was changed from 
10 kpc to 2 kpc. Thus, we only show the spectral fitting 
with Dbh = 10 kpc (see Figure |6|). The BH spin can 
take values in a wide range, from —0.998 to ~ 0.7. The 
unabsorbed flux in the 2-12 keV band are 2.7 x 10“® 
erg cm“^ s“^ and 2.9 x 10“® erg cm“^ s“^ for the non- 
flare and flare spectra, respectively. Using the average 
flux over the non-flare and flare stages, and assuming a 
typical Mbh of 10 Mq and Dbh = (2 — 10) kpc, the 
source luminosity is (1 — 34) x 10^® erg s“^ and the Ed¬ 


dington fraction {L/LEdd) is 0.1% - 2.7%. While the 
upper part of the L/LEdd range would not be unusual 
for a soft state, the lower part o f the range is low for a 
soft state (e.g., I Yu fc 120091 ). and this may favor a 
source distance closer to 10 kpc than 2 kpc. 

Previously, usin g the BeypoSAX M EGS spectra from 
the 2002 outburst. iMiller et al.l (1200911 measured the spin 
of the black hole in XTE J1908-I-094. The thermal emis¬ 
sion was not detected in these spectra; thus, they used 
the reflection component to constrain the spin and re¬ 
ported a value of 0.75 ± 0.09. If we fix the spin at 0.75 
and set the other parameters to be the same as for model 
5, the quality of the spectral fit is still good. Other than 
the BH mass being larger, the parameters are similar to 
those of model 5. However, given the large uncertainties 
in the spin, distance and inclination, it is impossible to 
constrain Mbh with our current data. The inclination 
measurement is independent of the BH spin that we as¬ 
sume with a value of ^ 30° — 40°, similar to i = 45° ± 8° 
reported bv IMiller et al.l (1200911 . 

Although these models containing the disk, the power- 
law (Comptonization) and the reflection components fit 
the NuSTAR spectra well, upon closer inspection, we 
find a small bump in the residuals near 8-9 keV (see 
Figure [S|). A similar feature is also obs erved in some 
other NuSTAR spectra, such as Cyg X-1 (|Tomsick et al.l 
1201411 . Adding a Gaussian emission line with Dcent 8.2 
keV and cr ~ 0.3 keV, the spectral fits are improved 
with Ay^ ~ 16, and the key parameters change only 
slightly. The line feature is likely related to a combi¬ 
nation of iron K/3 and nickel emission, nei ther of which 
are included in the reflionxJic model (|Walton et al.l 
in preparationll . 


4. DISCUSSION 

We have presented NuSTAR and Swift observations of 
XTE J1908-I-094 during its 2013 outburst. Similar to the 
two previous outbursts, the source was first detected in 
the hard state, then went into the soft state and decayed 
rapidly afterwards. The time interval between the 2013 
outburst and the last one is about 10 years, which is much 
longer than that between the two previous outbursts. 

The NuSTAR light curves show a flare of ~ 40 ks dura¬ 
tion. Fitting the spectra with the two components com¬ 
bining model of diskbb plus pegpwrlw prior to, during 
and after the flare, we found the power-law component, 
rather than the disk component, exhibited major changes 
during the flare. The power-law was softer and brighter 
during the flare but seems to be stable in the stages prior 
to and after the flare. A possible scenario is that there 
was an injection of high-energy particles (perhaps due 
to a jet ejection or shocks in the accretion disk) dur¬ 
ing the flare; thus, the flux of the power-law component 
increased, and the power-law index varied. If we keep 
the power-law model constant during the whole obser¬ 
vation and add another power-law model in fitting the 
flare spectrum, the extra emission is found to have a 
photon index of F = 2.61 ± 0.04 with a 3-79 keV flux of 
^ 1.8 X 10“^® erg cm“^ s“^. 

Jet ejections are not unusual in Galactic X-ray bi- 
naries. Other sourc e s, su c h as CHS 1915-1-1 05 (e.g. 
iMirabel fc RodriguezI 11994 : iFuchs et al.l I2003D . GRO 


J1655-40 (e.g . Tingav et al. I1995D and Cwnus X-1 (e.g. 
iStirling et al.ll200lD . also show ejection events. There are 

















































Tao et al. 


at least two types of ejections: one is the discrete out¬ 
flow usua lly appearing in the hard-to-soft sta te transi¬ 
tion fe.g.. lFender et al.1l2004lCorbel et al.ll2004ll . and the 
other is the compact jet occurring in the hard state and 
disappearing in the hard-to-soft state tran sition (e.g., 
iFender et al.lll99^ iCorbel et alJ[200nl 1200, Hi . The radio 
flux of XTE J1908-I-094 showed a signific ant increase be¬ 
tween 2013 November 5 and November 6 (|Rushton et alJ 
l2013b[l , which was 2-3 days be fore the NuSTAR obser¬ 
vations. iRushton et al.l (|2013bl l and iCoriat et al.l (|2013f ) 
suggested that the source ejected some optically thin 
radio-emitting plasma during the period. In fact, the ra¬ 
dio flux peaked during the NuSTAR observation. Also, 
the radio polarization measurements and the radio spec¬ 
trum are consistent with the discrete ejection interpreta¬ 
tion (jCnrran et al.l[2?)T^ . Moreover, we note that radio 
flares are a ccompanied by X-r ay flares in some X-ray bi¬ 
naries (e.g. IWilms et al.ll2007h . Based on the facts that 
the source was in a transition from the hard state to the 
soft state and the radio flare was apparent, such X-ray 
ejections would not be surprising. The X-ray flare ob¬ 
served by NuSTAR may have been caused by a discrete 
ejection. 

Although a couple of faint X-ray jets were detected 
from microquasars a few years after the ejection (e.g., 
ICorbel et al.l 1200^ 120051) . the X-ray emission produced 
by the plasmoid ejection may not be enough to explain 
the observed flux of XTE J1908-I-094. Thus, sudden 
changes of the temperature or the size of the corona 
may be anothe r scenario. We used a Comptonization 
model comptt (|Titarchukl I1994D to replace the power- 
law model in model 1, and untied different combina¬ 
tions of parameters in the non-flare and flare spectra. 
If the plasma temperature, fcTe, and the 3-79 keV flux 
in units of 10“^^ erg cm“^ s“^, A^comp, are allowed to 
be different for the different spectra, we obtain kT^ = 
500^140 keV and iVcomp = 377 ± 4 from the non-flare 
spectra, and kT^ = 306^77 keV and N^omp = 526 ± 5 
during the flare, with x^/dof = 2660.5/2098; if, instead 
of the temperature, the plasma optical depth parame¬ 
ter, r, is allowed to be free, we obtain t = 0.030to oo 2 
and A/;omp = 376 ± 4 from the non-flare spectra, and 
r = 0.010^0°^^ and A^comp = 527 ± 5 during the flare, 
with x^/dof = 2658.7/2098. We note that the ejection 
might remove material and cause a drop in the optical 
dept h. If the corona is act ually part of the jet, such as its 
base (jMarkoff et al.l[^05t) . the two explanations that we 
discuss (an ejection or a change in the coronal properties) 
might be related. 

During the hard state of the 2002 outburst, a broad line 
feature with an average energy oi E = 5.73 ± 0.09 keV 
and a line w i dth of a = 1.11 ± 0.31 keV was observed by 
iGogfis et al.l (j2004D . The feature disappeared in the soft 
state, whereas it reappeared when the source later en¬ 
tered into the hard state. The flux of the line component 
is str ongly linked to tha t of the power-law component; 
thus, IGogfis et al.l (|2004D suggested that the line feature 
might be the Fe Ka line from the reprocessing of the hard 
X-ray photons by cooler material close to the central ob¬ 
ject. In order to search for the line feature o ver the whole 
2013 o utburst, we used the same model as IGogfis et al.l 
(|2004D to fit the Swift spectra. However, for most obser¬ 


vations, the line feature is not remarkable, and the two 
component model containing the diskbb and power-law 
components could also fit the spectra successfully (Fig¬ 
ure [21 d). This may be due to the lower throughput of 
Swift/XRT above 6 keV, making the line feature unde¬ 
tectable. Moreover, in several observations of the soft 
state, the spectra show a possible Fe Ka line feature, 
and this is further confirmed by the NuSTAR observa¬ 
tion (see Figure |4] and Figure [5] where there is an iron 
line and hard X-ray bump). The Gaussian emission line 
models could fit the line feature, with model 2b (free 
A^h) providing a better fit (see Figure [5] and Table [2]). If 
a moderate A^h of 2.5 x 10^^ cm“^ or 4.3 x 10^^ cm“^ 
is used, Ecent agrees with the energy range of iron emis¬ 
sion; if A^h is allowed to be free, Ecent is well below this 
energy range, in which case the emission line may be red- 
shifted due to the gravitational effect. Regardless of the 
value of A^h, the line widths are about 1-2 keV , which 
are similar to those rep orted by lin’t Zand et al.l l)2002D 
and IGogfis et al.l (|2004l) . 

The reflection component can also be well fitted by 
the reflionxJic model, although the folding energy is 
not well constrained. As would be expected due to the 
stronger power-law flux during the flare, the ionization 
parameter, during the flare is larger than for the non¬ 
flare spectra. 

Adding a relativistic blurring model (model 4), 
relconv, provides only a small improvement on the qual¬ 
ity of the fit to the spectrum, and the parameters of 
model 4 and model 5 (replacing the multi-temperature 
disk-blackbody by the kerrbb model), agree with those 
of model 3. The reflection covering fractions, calculated 
from the ratio in 20-40 keV flux between the reflection 
and the power-law component, are 1.1-2.2 for the non¬ 
flare spectra and 1.6-3.9 for the flare spectra based on 
the different models. We note that the covering frac¬ 
tions are larger than 1, which indicates that the X-ray 
emission may come from closer to the black hole and the 
relativistic effec ts are stronger so that the light is gravi- 
tationally bent (iMiniutti &: Fabiaiil 12004 iMiniutti et al.l 
I2004D . Based on the relconv parameters, all possible 
values for the spin of the black hole (—0.998 to 0.998) 
are allowed when the inner radius is fixed to the ISCO, 
which means that the spectra may be extremely blurred, 
with a maximal BH spin, or somewhat less blurred with a 
retrogr ade disk or wi t h the inner disk being ionized. Fol¬ 
lowing iDauser et al.l (|2014D . the reflection fraction can 
give a lower limit on the black hole spin when assum¬ 
ing a lamppost geometry (i.e., a point-like corona above 
the spin axis of the BH); in that case a covering fraction 
above 1.6 implies a spin greater than 0.6. 

5. SUMMARY AND CONCLUSIONS 

NuSTAR and Swift observed XTE J1908-I-094 during 
its 2013 outburst. The Swift monitoring observations 
show that the source reached the soft state very close to 
the time that NuSTAR observed the source. A flare with 
a duration of ~ 40 ks appears in the NuSTAR light curve, 
peaking at ^^40% above the non-flare level. When fit¬ 
ting the non-flare and flare spectra with two-component 
models, consisting of diskbb plus pegpwrlw or diskbb 
plus comptt, we found that the power-law component 
(or the Comptonization component), rather than the disk 
component, went through great changes during the flare. 
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Changes of the corona, including variations of its temper¬ 
ature or its size, or the ejection of hot plasma, are two 
possible and potentially related scenarios for the flare. A 
broad iron line feature with cr = 1 — 2 keV is observed 
in the NuSTAR spectrum, which motivates a spectral 
model that combines a thermal disk, a power-law and a 
reflection component, providing a good fit to the spec¬ 
trum. Although the broad iron line provides evidence 
for relativistic blurring of the reflection component, we 
are not able to constrain the BH spin in the spectral 
fits, and all possible spin values, from —0.998 to 0.998, 
are allowed. The strong reflection component requires a 
covering fraction in excess of 1.0, which may be explained 
if light bending by the BHs gravitational field enhances 
the flux incident on the inner disk and suggests a spin 
larger than 0.6. 
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